Conductivity of electronic liquid-crystalline mesophases 
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We investigate the connection between the transport properties and the thermodynamics of elec- 
tronic systems with a tendency to form broken-symmetry mesophases evocative of the physics of 
liquid crystals. Through a hydrodynamic approach to the electronic transport in inhomogeneous 
systems, we develop a perturbative expansion for the macroscopic conductivity to study the trans- 
port of two-dimensional smectic and nematic phases. At the fluctuation induced flrst order phase 
transition expected for the smectic to isotropic transition, a jump in the macroscopic conductivity 
is predicted, with a directional dependence that reflects the fluctuation spectrum of the order pa- 
rameter. When elastic fluctuation modes melt the smectic phase into a nematic phase, the resultant 
nematic order parameter is shown to be linearly proportional to the conductivity anisotropy. We 
also outline qualitative comparisons with recent experimental works on strongly correlated materials 
that show evidences of electronic liquid-crystalline mesophases. 
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I. INTRODUCTION 

Recent investigations on strongly correlated electron 
systems suggest the emergence of inhomogeneous charge- 
ordered phases reminiscent from the smectic and nematic 
states commonly found in Hquid crystals Such an 
analogy makes reference to the broken symmetries of 
each phase: while the nematic phase breaks only rota- 
tional symmetry, the smectic state also breaks the trans- 
lational invariance along one particular direction. Such 
electronic mesophases can be envisaged as "fluctuating" 
and static charge stripes, respectively. Numerous exper- 
imental works provide evidences for the occurrence of 
smectic charge-ordering in nickelates 0, manganites 
Q and LBCO cuprates Q, as well as indications of the 
existence of nematic order in quantum Hall systems [61, 
ruthenates 0, LSCO and YBCO cuprates 

The order of these mesophases is directly manifested 
in anisotropies of the electric transport as well as in the 
magnetic or charge response of the system, similar to the 
anisotropic optical and hydrodynamic properties of Hq- 
uid crystals, such as birefringence and viscosity An 
interesting question is the explicit connection between 
the order parameters of these inhomogeneous mesophases 
and their transport properties. While it is natural to ex- 
pect that the resistivity anisotropy is directly linked to an 
order parameter that breaks rotational or translational 
symmetry, it is less obvious how to connect transport 
properties to spatial or temporal fluctuations and corre- 
lations of the system. 

In this paper, we investigate the transport properties 
of inhomogeneous electronic mesophases in the hydrody- 
namic transport limit, where collisions dominate. In this 
regime, electron transport can be understood as a diffu- 
sive process with inhomogeneous, potentially time depen- 
dent diffusion constant D (x,i). We use a coarse grained 
description and assume that the dephasing length scale 
lo, beyond which diffusive transport sets in, is small com- 
pared to the length scale ^ on which the inhomogeneities 



of the electronic mesophases vary. The condition ?o ^ C 
is expected to be valid close to flnite temperature phase 
transitions. This is certainly true close to second order 
transitions where £, diverges. However, even if the tran- 
sition is weakly flrst order, we still expect ^ to be larger 
than the lattice constant a while becomes compara- 
ble to a due to strong inelastic scattering close to the 
transition. 

As will be discussed in detail below, our formalism 
for the conductivity of inhomogeneous mesophases has a 
close connection to the theory of random resistor network 
(RRN) Previously, the classical RRN was appHed 
to explain transport properties of composite fllms 
manganites [l^ and silver chalcogenides jT^. Moreover, 
numerical simulations of correlated versions of this net- 
work were carried out in the contexts of manganites [l7| . 
disordered electronic nematic phases in cuprates and 
flnite temperature Mott transitions [l^ ■ The present ap- 
proach is an analytical theory for transport in electronic 
mesophases valid for inhomogeneities that are small in 
amplitude but well correlated in space. It can also be 
considered as an analytic theory for resistor networks 
with correlated local conductivities [131 ■ Here, the spatial 
correlations between the microscopic resistors is directly 
connected to the correlations of the order parameter de- 
scribing the mesophase, and not to an arbitrary distribu- 
tion function j2l|. In the formalism that we develop, not 
only the order parameter mean value, but also its rele- 
vant fluctuations are explicit related to the macroscopic 
d.c. conductivity through a perturbative expansion. 

By applying this general hydrodynamic transport 
model to the static charge striped phase, we show that 
while the conductivity measured perpendicular to the 
stripes probes mainly the order parameter amplitude, 
the conductivity measured parallel to them is particu- 
larly sensitive to the fluctuation spectrum. Thus, from 
transport measurements, it is in principle possible to 
obtain information about the microscopic character of 
the anisotropic mesophase, such as the flrst moments of 
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the Boltzmann distribution function, for example. As 
a specific realization of the thermodynamics of the elec- 
tronic smectic state, we consider the Coulomb frustrated 
Ising model, first introduced by Emery and Kivelson 
in the context of doped Mott insulators that present 
high-temperature superconductivity [13, US]- We obtain 
an analytic expression for the d.c. conductivity by us- 
ing the self-consistent mean-field solution of the effec- 
tive Ginzbur g-La ndau Hamiltonian, called the Brazovskii 
Hamiltonian [2j|. We then show that the fluctuation in- 
duced flrst-order transition from the isotropic liquid to 
the smectic phase is manifested as an anisotropic jump 
of the conductivity, whose sign brings information about 
the fluctuation spectrum of the order parameter. Quali- 
tative comparisons to experimental data showing jumps 
in the resistivity of nickelates are also outlined. 



In our theory, the investigation of the transport proper- 
ties of the electronic nematics focuses on the role of ther- 
mally excited modes. Thus, it is complementary to the 
analysis of quantum modes as discussed in the recent lit- 
erature [i^. Based on the work of Toner and Nelson [2^, 
we also describe the high-temperature nematic phase as 
a smectic phase melted due to the elastic fluctuations of 
stripes. We show explicitly that there is a temperature 
range where the macroscopic conductivity anisotropy is 
linearly proportional to the nematic order parameter, as 
expected from symmetry considerations ^27)], and deter- 
mine the prefactor that connects transport properties 
and nematic order. The same linear relation, including 
the prefactor, was also observed in numerical simulations 
of a disordered electronic nematic phase at zero temper- 
ature , suggesting a close connection between the two 
rather different approaches. Finally, a very recent experi- 
ment regarding YBCO showed a remarkable resemblance 
between the spectral weight of the low-energy anisotropic 
spin fluctuations, obtained through neutron scattering, 
and the resistivity anisotropy, obtained from transport 
measurements II 11 1. 



The remainder of the paper is organized as follows: in 
section m we provide a detailed derivation of the hydro- 
dynamic model for the diffusion of an electron in an inho- 
mogeneous medium, which will be used throughout the 
work. Section Hnl is devoted to the application of this for- 
malism to the transport properties of an electronic smec- 
tic phase. Not only do we outline very general properties, 
but we also obtain speciflc results after describing the 
charge stripes thermodynamics by the Brazovskii model. 
In section llVt the d.c. conductivity of an electronic phase 
with nematic type of order is investigated through the hy- 
drodynamic transport model. Such a state is described as 
a smectic phase melted by thermally excited elastic fluc- 
tuation modes of the stripes. Comparisons to other ap- 
proaches as well as to recent experiments involving doped 
transition metal oxides are delineated. Finally, section Ivl 
is devoted to the flnal remarks and acknowledgments. 



II. DIFFUSIVE TRANSPORT IN 
INHOMOGENEOUS MESOPHASES 

Let us consider the diffusion of an electron in an arbi- 
trary inhomogeneous medium. We start from the conti- 
nuity equation 



dp (x, t) 
dt 



V-j(x,t) - 0. 



(1) 



connecting the charge and current densities. The macro- 
scopic conductivity of the medium relates the mean 
current density to the external electric fleld — VC/oxt (x,t) 



(ja (X, t)) 



<j^0 {t - t') V/3[/cxt (x,t') dt' , (2) 



yielding, in the Fourier space 

(p(k,w)) = -i—^ ^ 



C^cxt (k,w) . 



(3) 



Next, we assume that locally the diffusive relation 

ja (x, t) = -XoDap (x,i) V/3C/10C (x,i) (4) 

holds for the current density in terms of the local electri- 
cal fleld — V?7ioc (x,t), where Da^ (x,i) denotes the inho- 
mogeneous and possible time dependent local diffusion 
coefficient. Within linear response, the uniform charge 
susceptibility xo connects the local charge variations with 
the difference between the local and external potentials 



P (x, t) = Xo [U\oc (x,i) - J7oxt (x,i)] 



(5) 



Combining the last two expressions, we are able to re- 
late the external and local potentials via 



with inverse diffusion operator 
This yields, for the charge density 



V„ [D^p (x,<) Vp] . 



P (x, t) = Xo 



4-' 



C^oxt (X,t) 



(6) 



(7) 



(8) 



After taking the conflgurational average over p(x, i), 
we can Fourier transform its mean value and obtain, by 
comparing to Eq. [3] 



cra/3 (w) = lim TTir^o [~i(^G (k, uj) ~ 1] 



(9) 



where G (k, to) is the Fourier transform of the average of 
the differential operator ([7]). In the d.c. limit, this flnally 
yields 

<^ai3 = Xo lim lim ^-^Re [G (k, uj)] . (10) 
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Using the Einstein relation [28|, we can identify 
the anisotropic macroscopic diffusion coefficient as 
lim lim r^T^Re [G (k, a;)l. Therefore, it is clear that, 

even if the tensor of the local diffusion coefficient behaves 
as (x, i) = D {:x.,t) Sai3, the global diffusion coeffi- 
cient can be anisotropic, as long as D (x,i) is a function 
of the direction of x. 

This hydrodynamic formalism has a one to one anal- 
ogy with the theory of random resistor networks (RRN) 
[23|. In its most elementary form, one considers two resis- 
tors with conductivities a a and as, randomly distributed 
with probabilities p and 1—p over the links of a network. 
An external electric potential Ucxt is then applied in each 
site through a local capacitor with specific capacitance C. 
As shown by Stephen j2§|, the macroscopic conductivity 
is given by a result identical to EqfTOj with C playing the 
role of the charge susceptibility while the local conduc- 
tivities of the RRN correspond to XoD (x). In the RRN 
problem, the average is performed over the distribution 
function of the resistor network, which is frequently as- 
sumed to be a binomial distribution, characterized by the 
probability p. 

The key difference between our approach and the RRN 
theory is that the distribution function for the local diffu- 
sion coefficient is determined by the distribution function 
of the order parameter. Let the order of an inhomoge- 
neous electronic nematic or smectic state be character- 
ized by a scalar density field p (x, t). Here, p (x, t) is the 
deviation of the the coarse grained electron density from 
its mean value. We then assume a simple connection 



Dap (X,i) = Da[3 [p(x,t)] 



(11) 



between the spatially varying diffusion coefficient and the 
electron density of the electronic mesophase, meaning 
that the temporal and spatial variations of the diffusion 
coefficient are determined solely by those of p(x,t). In 
our cases of interest, the ordered inhomogeneous state is 
characterized by an order parameter that varies in space, 
alternating between p > and p < 0. For instance, in 
electronic smectics, p > (p < 0) denotes a hole rich 
(poor) coarse-grained region. Since each of these regions 
has its own conducting properties, we can associate dif- 
ferent local conductivities to each of them. 

For weakly inhomogeneous systems, we expand 
Da/3 (x,t) relative to the homogeneous state, where 
p (x, t) — 0. Hence, we propose the following specific 
form for Da/i (x,i) of Eq. lfTTI) 



Dap (x,<) = XOCTO [1 + gP (x,t)] Sap 



(12) 



where ctq is the microscopic conductivity mean value and 
g is the coupHng constant measuring the contrast be- 
tween the conductivities of distinct regions. The physi- 
cal meaning of these two parameters can be better vi- 
suaHzed considering the limit of an ordered homoge- 
neous phase where p = ±po. Denoting the conduc- 
tivity of the saturated homogeneous p > (p < 0) 



phase as cr> (cr<), we have = ((T> -I- cr< ) /2 and 
9 = Po^ ((t> — cr^) / (cr> -I- (T<). A relation similar to Eq. 
[121 was considered in the context of the Mott transition 
by Papanikolaou et al. [l^. However, in that case it was 
necessary to include a next order term cx p (x,t) in the 
expansion, EqfT2l originated from the contribution of in- 
terface scattering. Such effects are not included in the 
present work. 

Now, we particularize our analysis to the situation 
where the inhomogeneities of the system are small, 
such that (7^1. Thus, the diffusion differential op- 
erator (O can be perturbatively expanded as G = 



^5" (GoT>)"Go, with 



n=0 



V 



dt 

DoV ■ [p (x,t) V] , 



(13) 



where the bare diffusion coefficient is Dq = xo^'o- In 
order to obtain a perturbative expansion for the conduc- 
tivity through Eq. [lOl we need to perform an average 
over p before taking the Fourier transform of V. There 
are two opposite limits one can consider: in the first one, 
which we shall call quenched limit, the electron diffuses 
faster than the field fiuctuates, probing a frozen config- 
uration of the order parameter. Hence, after taking the 
proper Fourier transforms, we obtain, to second order in 
9 



l+.g(p(k = O,c^ = 0)) 



(14) 



d'^kduj- 



k'^ [ria ■ i^ 



k^) 



(p(k,w)p(-k, - uj)) 



where fia is the direction taken to measure the conduc- 
tivity and (• • •) denotes the proper average over the or- 
der parameter, which can be the usual thermodynamic 
average, for example. Similarly we could also consider 
transport in nonequilibrium configurations such as glassy 
states; then, (• • •) refers to the corresponding dynamic 
average of the nonequilibrium configuration under con- 
sideration [sol. [3l]|. For future reference, we rewrite the 
previous formula for the special situation where the field 
is static 



ao[l + .g(p(k = 0)) 



(15) 



'9 



d''k(ha--k\ (p(k)p(-k)) 



In the opposite Hmit, which we shall call annealed 
limit, we consider a field that fiuctuates much faster than 
the electron diffuses. Therefore, it is legitimate to replace 
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the order parameter by its mean value p (p) in the mi- 
croscopic conductivity expression lfT2|) . yielding 

(Tea = CTO [1 + .g (P (k = 0)) (16) 

-g2|d'^fc(n„.k)'(p(k))(p(-k)) . 

It is clear that the anisotropic character of the conduc- 
tivity is not manifested until second order in g. Moreover, 
this second order term reflects a fundamental difference 
between the quenched and annealed limits. While in the 
latter the conductivity behavior is dictated by the order 
parameter amplitude, in the former we note a tight con- 
nection between the macroscopic conductivity ctqq and 
the correlation function (p(k)p(— k)), resembling other 
similar relations found in condensed matter systems, such 
as between the scattering cross section and the thermo- 
dynamic correlation functions. 

III. TRANSPORT IN THE SMECTIC PHASE 

Using the formalism developed in the previous sec- 
tion, there are some very general statements we can 
make about the conductivity of an electronic smectic 
mesophase, independent on the speciflc model under con- 
sideration. Hereafter, we consider that the static order 
parameter p(x) describes local fluctuations of the charge 
density and that the electron difi^uses much faster than 
the fleld changes (quenched limit). Since the system 
is usually electrically neutral, such that {p (k = 0)) = 
y J (p(x)) d'^x = 0, the lowest non-vanishing correction 
to the uniform conductivity in perturbation theory is of 
second order in the contrast. After writing the correla- 
tion function as 

{p (k) p i-k)) = C (k) + {p (k)) {p (-k)) , (17) 

where C (k) is the connected correlation function, we can 
split the second order term of the macroscopic conduc- 
tivity (mi) in two parts, aaa = do [l - 5^ (o'Lq + ^aa)] : 

= J d'*A:(n„.k)'c(k) 

= / d'k (n„ -kfip (k)) {p (-k)) . (18) 

In the case of interest here, the correlation function 
depends only on the modulus of the momentum, i.e. 
C (k) = C{k). Therefore, a'^^ does not depend on the 
direction taken to measure the conductivity, since iia 
will be integrated out over all directions. Hence, cr^^ is 
isotropic and depends only on the order parameter fluc- 
tuation spectrum 

a'^^ (X / d^kC (k) = C (x, x) = (x)) - (p (x))^ . (19) 



Meanwhile, the term a'^^ is proportional to the aver- 
aged order parameter. For small amplitude static stripes, 
it is reasonable to assume that the order parameter mean 
value is strongly anisotropic and has a pronounced peak 
along the modulation direction qo, such that (p(k)) ~ 
A(5 (k — qo). Higher harmonics (entering as (5 (k — 2qo) 
or (5 (k — 3qo) etc.) only matter once the amplitude of 
the inhomogeneity becomes large, i.e. p (x,i) becomes 
large compared to the mean electron density. It follows 
that, unHke a'^^, the term is anisotropic and depends 
on the relative angle 6a between the external electric field 
and the modulation vector through 

cx cos2 Oa . (20) 

Thus, if the conductivity is measured perpendicular to 
the stripes (n^ |j qo), the order parameter and its fluctu- 
ation spectrum will be probed. However, if the current 
is applied parallel to the stripes (ria -L qo), the measure- 
ment will be sensitive solely to the fluctuations of the or- 
der parameter. Clearly, the transport quantity that only 
probes the order parameter amplitude is the anisotropic 
conductivity a± — . 

Another interesting situation to investigate is when the 
system undergoes a phase transition to the state with no 
broken symmetries, the isotropic liquid. For instance, let 
us consider a flrst order transition from the liquid to the 
smectic phase: if the conductivity is measured parallel to 
the stripes direction, it is expected that cr'aa vanishes, and 
the isotropic part tr^^ will dictate the behavior of the d.c. 
conductivity. Since fluctuations of the order parameter 
are usually larger in the liquid disordered state, we expect 
from Eq. [I9]that, close to the transition, the conductivity 
of the liquid phase will be smaller than the one of the 
ordered smectic phase, a^^"^ < a^™'^. 

Instead, if the conductivity is measured perpendicular 
to the stripes direction, the contribution of the term tr^^ 
will be maximum. Unlike its counterpart a'^a, this term 
is greater in the ordered side, since the order parameter 
amplitude vanishes in the liquid phase. Hence, the ques- 
tion of whether a'^'^ or (t^™'^ is larger close to the phase 
transition depends on the ratio between the jumps of the 
fluctuations and of the amplitude of the order parameter. 
For strong flrst order transitions, we expect the latter to 
be more significant, meaning that > a'^'^. 

Let us now apply the hydrodynamic transport formal- 
ism to a specific model for the thermodynamics of the 
electronic smectics. In particular, we focus on doped lay- 
ered transition metal oxides, which are believed to be 
well described as 2D doped Mott insulators [s^l- In these 
compounds, the stripes appear as a compromise between 
a short-range tendency of the doped charges to phase 
separate from the antiferromagnetic background spins 
and the long-range electrostatic repulsion between alike 
charges [121 ■ A simple model that attempts to capture 
these properties is the Coulomb-frustrated Ising model 
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27r 



(21) 



where p(x) > denotes a hole-rich region (larger con- 
ductivity) while p(x) < denotes a hole-poor region 
(smaller conductivity), with / {p{x)) d'^x — 0. As usual, 
To denotes the reduced temperature and u is an effective 
parameter. While the first term of the above Ginzburg- 
Landau Hamiltonian expresses the tendency of phase sep- 
aration, the second term accounts for the frustration in- 
troduced by the Coulomb interaction. This Hamiltonian 
is minimized in the Fourier space for a non-zero wave 
vector whose modulus is go- Thus, considering the low- 
energy fluctuation modes, we can expand Eq. [21] around 
its minimum to obtain the effective functional: 

n ^ ^1 d^xd^x'p{^)c-'{^,^')p{^) + l I d^xpH^) 

pik-(x-x') 

d^k . (22) 



Co-i(x,x') 



To + {k - qo) 



Note that although tq in Eq. [22] is shifted with respect 
to To in Eq. [211 we keep the same symbol to simplify the 
notation. This model was first studied by Brazovskii in 
the context of cholesteric liquid crystals [l^l, and has 
been employed to describe a variety of other physical 
systems with inhomogeneous states, from "hard matter" 
ones, such as dipolar ferromagnets [33jj3J|, to "soft mat- 
ter" s yste ms like diblock copolymers [3q| and microemul- 
sions f36|. The striped phase with modulation vector qo, 
which has smectic order, is a thermodynamic stable state 
of the model. In two dimensions, extra external poten- 
tials are necessary in order to ensure the stability of this 
smectic phase against elastic fluctuations - like crystalline 
flelds [37| or pinning centers, for example. 

Since we are interested in the main qualitative trans- 
port properties of the system, we use the self-consistent 
mean-fleld solution of Eq. [22] The details of this method 
are presented elsewhere (24 Isl. Issl . l39| . It predicts a 
fluctuation induced flrst order transition from the smec- 
tic state (p(x)) = 2A cos (qo • x) to the liquid disordered 
phase. The temperature for which their free energies be- 

/ r- n2/3 

come equal is Tq w — (3v2(?ouj , which we will con- 
sider the transition temperature hereafter (this is valid 
for an adiabatic change of the temperature). At this 
point, there is a jump in the correlation length of the 
system ^, as sketched in flgure[TJ We consider the stripes 
frozen in the diffusive electron reference frame; hence, 
substituting the self-consistent solution of l(22|) in Eq. [15] 
yields, for the conductivity of the liquid disordered and 
smectic phases 



liq 



(23) 



2.0 



1.5 



1.0 



Figure 1: Correlation length ^ (in units of the correlation 
length 5* of the smectic phase at the reduced transition tem- 
perature Tq) as function of the reduced temperature tq. 




Figure 2: Conductivity a± measured perpendicular {Oa = 0) 
to the stripes direction (solid line) as function of the reduced 
temperature to (in units of the transition temperature modu- 
lus |to I). An intrinsic insulating behavior was considered for 
(70 (dashed line). 



(Tq 



cos 



where ^ is the correlation length, ^* is its value at the 
transition temperature for the smectic phase and 6a is 
the angle between the direction taken to measure the con- 
ductivity and the direction of the stripes modulation, qo. 
As expected, the conductivity is intrinsically anisotropic 
and it is greater when measured parallel to the stripes 
{Oa = f ) than when it is measured perpendicular to them 
{6a = 0), cr^™'^ > CTjJ"'^. Moreover, this anisotropy in- 
creases as the temperature is lowered due to the decrease 
in Another manifestation of the thermodynamic be- 
havior in the transport properties is the occurrence of a 
jump in the conductivity at the transition temperature, 
as shown in flgure[2j 

Not only does the module of the jump Act {6a) — 
<Taa I * ~ '^aa'^lr'- depend On the angle between the exter- 
nal current and the stripes modulation, but also its sign. 
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Figure 3: Conductivity jump Act (9) as function of the angle 
between the modulation direction qo and the direction taken 
to measure the conductivity (in degrees). 



due to the coupling to the lattice, in accordance with 
neutron-diffraction measurements that verified a decrease 
in qo ji^l- A possible explanation is given in the con- 
text of the present model: considering a reduction in the 
wave vector modulus q'^ = rjqo, rj < 1, and solving the 
resulting self-consistent Brazovskii equations, we obtain 
fjsmc ^g/^-j ^ CT^™^ (gp) for T] close to 1 and for any temper- 
ature below the transition one. Although other mecha- 
nisms could also explain this decrease in the conductivity, 
like a temperature dependent contrast g, for example, it 
is worth pointing that our model is consistent with this 
experimental observation. 

IV. TRANSPORT IN THE NEMATIC PHASE 



In figure m we plot Act {9a) as a function of 9a and see 
that while its sign is positive when the conductivity is 
measured perpendicular to the stripes, it becomes nega- 
tive when the measurement is made along the stripes. 

This is in agreement with the general discussion carried 
out in the beginning of the present section. The negative 
sign is a consequence of the fact that the conductivity 
measured parallel to the stripes {9a = f ) probes the 
fluctuations of the order parameter, which are greater in 
the disordered phase, as it can be noted directly from the 
behavior of the correlation length presented in figure [H 
When the conductivity is measured perpendicular to the 
stripes {9a — 0), even though fluctuations still contribute 
to the conductivity, there is an extra contribution coming 
from the order parameter amplitude, which is non-zero 
only in the ordered phase. The positive sign is a balance 
between these two contributions, meaning that the latter, 
proportional to is greater than the flrst, proportional 
to 

In the LTO phase of the layered nickelates 
Ndi_xSrxNi04, a sudden upturn in the resistivity was 
observed at the charge-ordering temperature 0, 01 • The 
external current in the experiment of Hiicker et al. Q 
was appHed perpendicular to the direction of the stripes 
and this sudden change could be the result of a jump 
rounded by disorder [40|. It would be interesting to re- 
alize an experiment where the conductivity is measured 
along the stripes direction in order to check if an op- 
posite change in the resistivity takes place. This would 
provide additional criteria to check the validity range of 
this hydrodynamic transport model and send Hght on the 
fluctuation spectrum in a state of mesoscale ordering. 

In another class of nickelates, the La2Ni04+5, experi- 
ments probing resistivity fluctuations indicated no jump, 
but also suggested a strong coupling to sample disorder, 
giving rise to slow stripe dynamics j4l|. The current 
was applied perpendicular to the stripes direction. For 
slow sample cooling, however, jumps were observed be- 
low the transition temperature, inside the smectic phase, 
with an increase in the resistivity. This was attributed 
to changes in the modulation wave vector modulus go 



In the previous section, we implicitly assumed the 
underlying presence of crystalline fields and/or pinning 
disorder in the system, in order to keep the smectic 
phase stable against low-energy fluctuation modes of the 
stripes. In the absence of such stabiHzing mechanisms, it 
is well known that, at any finite temperature, the clean 
two-dimensional smectic phase is unstable towards the 
formation of a nematic phase due to elastic fluctuations 
of the stripes walls \2^. Therefore, the nematic phase 
can be conceived as a melted smectic phase, where the 
stripes have no longer true positional order. Since two 
is the lower critical dimension of the resulting system, 
the 2D nematic state has only quasi-long range orienta- 
tional order, which is lost above the Kosterlitz-Thouless 
transition temperature Tkt where pairs of disclinations 
unbind, driving the system to the isotropic liquid phase. 

Fluctuation modes of stripe walls other than thermally 
excited ones can also melt the smectic phase. For exam- 
le, random-fleld disorder [l^ and quantum fluctuations 
]J are able to drive the system to a nematic state even 
at zero temperature. Alternatively, the quantum elec- 
tronic nematics can be envisaged not as a result of the 
quantum melting of the smectic phase, but as a conse- 
quence of quadrupolar Pomeranchuk instabilities of the 
Fermi surface of the isotropic liquid phase j2l,li3|- Here, 
we will not consider the role of quantum fluctuations nei- 
ther the presence of disorder, but will follow Toner and 
Nelson's description of the nematic phase, taking into 
account only thermally excited elastic fluctuations [i^ . 
This is fully consistent with our hydrodynamic approach 
where the dephasing length is assumed short as the sys- 
tem is considered at finite T. 

Without loss of generality, we consider stripes modu- 
lated along the x direction and substitute the order pa- 
rameter p(x) = 2Acos [go {x + u{x,y))] in Eq. [221 where 
the displacement field u(x) describes the elastic fiuctua- 
tions. The resulting elastic Hamiltonian is given, in the 
harmonic approximation, by 



"^clastic 2 



d^X 



du\' 
dx J 



+ A2 



(24) 
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with B = A^q^ and = 250- In two dimensions, ex- 
plicit calculations of the mean value and of the correlation 
function of the order parameter yield (p) = and 



(p(x)p(O)) (X 



exp 



9o 
B 



47tA 



1/2- 



for Ax ^ 



(25) 

meaning that the system has not true positional order. 
After taking into account the roles of two kinds of topo- 
logical defects, dislocations and disclinations, it follows 
that the elastic Hamiltonian describing the melted state 
is the Frank Hamiltonian 12611 



Wprank =\ \ ^ {^1 ' + [N X (V X N)]'} , 

(26) 

which describes the elastic properties of a 2D nematic 
liquid crystal. In the previous expression, the unit vector 
N denotes the nematic director and the Frank constants 
K\ and Kj, are functions of the elastic parameters B and 
A, as well as of the energy necessary to excite an isolated 
dislocation, Ed. 

The role of the topological defects is fundamental to 
the understanding of the thermodynamics of the result- 
ing nematic phase. As explained above, the unbinding 
of disclination pairs induces a Kosterlitz-Thouless tran- 
sition to the isotropic liquid phase. Moreover, since iso- 
lated dislocations have a finite excitation energy E^, a 
new length scale is introduced in the system below 
Tkt- Such a length scale denotes the correlation length 
of isolated dislocations that proliferate in the system as 
the temperature is increased, and is given by = n-^^"^ , 
where ud = exp (^—-^p^ is the dislocations density 

and d oc q^^ is the stripes mean width. For lengths 
greater than the system has the same properties of 
a nematic liquid crystal, while for smaller lengths there 
is smectic order decorrelated by elastic fiuctuations only. 
Therefore, one can consider the nematic phase as a set of 
finite size smectic blobs [1^, where each blob has orien- 
tational order given by the average of the orientation of 
the stripes enclosed, as sketched in figured 

We can now apply the formulae of the previous section 
to calculate the conductivity of an individual blob cr^^^. 
Inside it, although density fiuctuations remain small as 
long as To < Tq, elastic fiuctuations are relevant. While 
a diffusive electron takes the typical time tei ~ '^o^S.d 
to cross the blob, phonon fiuctuations propagate with a 
typical time scale tph ~ v~^^d, where w is a sound veloc- 
ity. Therefore, we have to substitute the order parame- 
ter mean value {p (k)) = A [5 (k - qo) -f (5 (k -h qo)] e"^, 

where W = 6^*=^^ ^ |^2i7r^A^^ is the Debye-Waller fac- 
tor of the finite-size blob, in the annealed limit of the 
macroscopic conductivity, Eq. [161 We obtain 




Figure 4: Schematic representation of the nematic phase orig- 
inated from the melting of the smectic state due to thermally 
excited elastic fluctuations. The blobs comprehend regions 
free of topological defects enclosing stripes that are decorre- 
lated by elastic fluctuations only. The director of each blob is 
indicated by a double arrow. 



-2^C0S2( 



(27) 



where 9a is the angle between the direction taken to 
measure the conductivity and the blob director (we are 
following the convention that the director is perpendic- 
ular to the stripes orientation inside each blob). This 
procedure is valid for intermediate temperatures only, 
where W is small and the elastic fiuctuations \5pei\ = 
\/e^^^^^ (pei) are not so large. 

To obtain the system's macroscopic conductivity, the 
fiuctuations of the nematic director inside each blob have 
to be taken into account. We can consider them slow 
in the electron reference frame. Thus, while taking the 
annealed limit to calculate the conductivity inside each 
blob, we have to take the quenched limit to calculate 
the conductivity of the whole system, averaging over all 
blobs. We obtain 



XX 

nem 

yy 



0-0(1-9 (cos 20) 
ao(l-Hg(cos2C)) 



(28) 



for the two opposite situations where the direction taken 
to measure the conductivity is parallel {{0a) = 0, x axis) 
or perpendicular {{9a) = 7r/2, y axis) to the mean ori- 
entation of the nematic director. In the previous expres- 
sion, ( — 9a — {9a) denotes the director deviation from the 
mean value, ao = ctq (l - g'^A'^e''^^) and g = g'^A'^e''^^ 
(to second order in the original g). 

It is remarkable that the nematic order parameter 
(cos2C) appears spontaneously in the conductivity ex- 
pression, even though it was not assumed any coupling 
between the microscopic conductivity and the nematic 
order parameter. Such a relation can be cast in a more 
formal way by considering the tensorial nematic order pa- 
rameter in two dimensions, {Sap) — S (^NaN/^ — ^Sap), 
with S = (2cos2C- !)• It is clear, from Eq. [28l that 
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the anisotropic conductivity can be written as cr"^™ ~ 

o-Q (1 -5(S'a/3))• 
Another direct consequence of Eq. [28] is that the con- 
ductivity anisotropy is proportional to the nematic or- 
der parameter, ct^^™ - crJI" = 2cjog^A^e-^^ (cos2C), in 
agreement with symmetry arguments [23|- By comparing 
this expression with the one referring to the conductivity 
anisotropy of the smectic phase, Eq. [231 we note that 
in the present case the anisotropy is much smaller due 
to the Debye- Waller factor e~^^. Moreover, this result 
agrees exactly with the relation 

(^) (tr^) H"-c> ... 

found for T = by Carlson et al. fl8?|, with the same pref- 

actors (in our case, r = ^'^[^^cos 2C)-^-ij ~ T^-'' "^"^ ^^^^ 
work, the disordered electronic nematics was mapped on 
the Random Field Ising Model (RFIM), and the conduc- 
tivity was obtained numerically through a random re- 
sistor network approach. Although the two systems are 
very different, since ours is clean whereas the one investi- 
gated by Carlson et al. is disordered, this alike outcome 
seems to be a consequence of the existence of clusters 
with short-range smectic order inside the nematic phase 
- the thermally excited blobs in our model and the dis- 
order generated nematic patches in the other work. 

Such a linear relation between a macroscopic quantity 
and the thermodynamic order parameter can be very use- 
ful in the investigation of the orientational ordering of 
nematic phases in the cuprates, specially through trans- 
port measurements. In fact, a recent experiment jTlj 
on YBa2Cu3 6.45 suggested a linear relation between 
the compound conductivity anisotropy and the spectral 
weight of the low-energy anisotropic spin fluctuations. 
Still regarding possible experimental implications, we 
also point that due to the large elastic fluctuations in- 
trinsic to the system's dimensionality, deviations from 
the averaged conductivity are expected for small systems, 
what would be manifested as noise in time series measure- 
ments. 

Before flnishing this section, an important remark re- 
garding the mean value (cos 2^) should be made. Strictly 
speaking, in two dimensions this average is zero in the 
thermodynamic limit, due to Mermin- Wagner theorem 
(see, for instance, Chaikin and Lubensky [l^l). How- 
ever, in real layered systems, flnite-size effects as well 
as small interlayer couplings are able to stabilize the ne- 
matic phase, granting a non-vanishing value for the order 
parameter. 



V. CONCLUDING REMARKS 

We presented a hydrodynamic transport theory that 
can provide important tools in the investigation of elec- 
tronic phases with smectic and nematic symmetries, ren- 
dering explicit relations between macroscopic transport 
quantities and the microscopic order parameter. Partic- 
ularly, we considered doped layered transition metal ox- 
ides, using the Brazovskii model to describe the thermo- 
dynamics of the low-energy charge modes of their smec- 
tic phases. The directional dependence of the sign of the 
conductivity jump was shown to be a manifestation of the 
Brazovskii fluctuation spectrum and a general character- 
istic of the hydrodynamic transport model, constituting 
an interesting criterion to decide on the applicability of 
the model to describe the observed static charge striped 
phases. 

The finite temperature electronic nematics was con- 
ceived as a smectic Brazovskii phase melted by thermally 
excited elastic fluctuations of the stripes walls, following 
the approach of Toner and Nelson [2a|. Not only does 
the nematic order parameter appear explicitly and spon- 
taneously in our formalism, but it is also shown to be 
linearly proportional to the conductivity anisotropy, fol- 
lowing the same relation found numerically by Carlson et 
al. [l8j for the case of the disordered electronic nematics 
at zero temperature. 

In the context of nematic phases, it would be interest- 
ing to investigate other excitations that are also able to 
melt the smectic phase in a nematic state, such as quan- 
tum and disorder induced fluctuations P, Moreover, 
additional studies on out of equilibrium properties of the 
electronic nematics, particularly on the connection be- 
tween fluctuations of the smectic clusters and transport 
time-series measurements, would also be desirable to pro- 
vide a richer picture of the problem. Finally, applications 
of the general formalism to systems other than doped 
transition metal oxides can also be envisaged, for exam- 
ple, for the low-temperature striped phases that appear 
in quantum Hall systems jH, [13, 53] and for the inhomo- 
geneous phases found in spin glasses and Mott insulators 
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